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The rotary inverted pendulum (RIP) has been used in various control 
application areas. This system can be represented as two degree of freedom 
(2-DOF), consisting of a rotating arm and rotating pendulum rod. RIP is an 
excellent example of designing a single-input multi-output (SIMO) system. 
Due to unstable RIP system dynamics, and its nonlinear model, multiple 
control techniques have been used to control this system. This paper uses 
integer and fractional order proportional integral-proportional derivative 
(PI-PD) controllers to stabilize the pendulum in the vertical direction. 
Constrained optimization approaches, such as the grey wolf optimization 
(GWO) methodology, are utilized to estimate the parametric values of the 
controllers. The simulation results showed that the fractional order PI-PD 
controller outperforms the integer order PI-PD controller with and without 


RIP system disturbance signal existence. A multiple results comparison has illustrated 
the superiority of fractional order controller over a previous work. 
This is an open access article under the CC BY-SA license. 
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1. INTRODUCTION 

Rotary inverted pendulum (RIP) is a good example of control theory verification in control 
engineering. It is an excellent model for controlling space booster rocket attitude, satellite attitude control, 
autonomous plane landing systems, airplane stabilization in unsteady airflow, ship cabin stability, segway, 
and humanoid robots [1]. Also, this system has new applications, such as energy harvesting systems, which 
are considered an effective solution to excerption the kinetic energy from rotary systems where a pendulum 
system is connected to an electromagnetic generator [2]. Inverted pendulum systems are underactuated, 
unstable, nonlinear, open-loop, and challenging to handle. Due to its inherent nonlinearity, it is an interesting 
subject from a control standpoint. However, establishing an accurate mathematical model of the process 
using differential equations is typically difficult [3]. 

In the control engineering field, there are various types of inverted pendulums. Directional inverted 
pendulums [4], [5], double inverted pendulums [6], and RIP [7] are the most common. This study investigated 
the RIP type, which designed as a single-input multi-output (SIMO) system [8]. The stability of the RIP in the 
upright direction steadies the focus of most rotary inverted pendulum research. Various controllers have been 
utilized to regulate the RIP. We can summarize several researchers’ work as: Sirisha and Junghare [9], used 
linear quadratic regulator (LQR), proportional integral derivative (PID), fuzzy logic, and Hoo controllers to 
regulate the RIP performance and perform desirable stability. They did not assure the robustness of the 
controller, and the pendulum angle took a long time to settle. Ali and Naji in 2016 [10] suggested a state 
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feedback controller discipline the RIP behavior. The rod and arm angle response took a long time to reach the 
desired output with and without the parameter variation existence, considering the arm set point to be 
achieved at a very small magnitude (5 degrees). Sunil and Manju [11] have been introduced fuzzy logic 
controller to stabilize the RIP in the vertical direction. However, they did not check the system response after 
exerting an external impact, also, the arm response has not been listed. Lanjewar et al. [12] applied a 
fractional order PID controller to perform satisfied stability for the RIP system. The disturbance rejection was 
unacceptable when they applied 1.2 Nm as an external force. Öksüz et al. [13] proposed a full-state feedback 
controller for the system using MATLAB for a linear dynamic model. The results show that the system states 
response speed, but no robustness test is applied to the system. 

The works listed above showd various control techniques utilized to stabilize RIP. The contribution of 
this work is summarizing the design of two closed-loop fractional-order proportional integral (PI) fractional 
order proportional derivative (PD) (FOPI-FOPD) control systems to control the RIP at the same time as a SIMO 
design approach, which is considered as one of the complex control systems. Rather than the complexity of this 
design, an optimization technique is accomplished to find the appropriate control system parameters. 
The modelling of RIP can be presented in this work using state-space representation, FOPI-FOPD, and integer 
PI-PD controllers that can be tuned by grey wolf optimization (GWO) technique used to control the RIP. The 
fractional order controllers (FOC) with appropriate performance outperform integer order controllers [14], [15]. 
The FOC like FOPD, FOPI, fractional order proportional integral derivative (FOPID), and FOPI-FOPD might 
perform better than integer-order [16], [17]. The quanser RIP is used in this work, as shown in Figure 1. 
An angle of a displaces the pendulum. However, @ is the rotating angle of the arm around the vertical axis. 


»@ ‘ 
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Figure 1. SRV-2 RIP model 


2. MATHEMATICAL MODEL OF RIP 

The rotating position of the RIP arm is obtained in Figure 1. The rod can be described as a lumped 
mass at mid length of the rod. The rod is moved by a angle. However, the 9 depicts arm displacement at the 
x-axis. Therefore, a mathematical model can be created by inspection at the rotational velocity of the 
pendulum rod center of mass (COM). The terms used in system model derivations are listed in Table 1 [18]. 


Table 1. RIP system parameters 


Variable Description 
L Rod to COM displacement (m) 
m Arm mass (Kg) 
r Arm length (m) 
0 Arm angle (radians) 
a Rod deviation (radians) 
h Rod COM from the ground distance (m) 


dex: Rod inertia about its COM (Kg.m’) 
Ve Rod velocity COM at the x-direction (m\s) 
V, Rod velocity COM at the y-direction (m\s) 


The velocity of the rod and arm can be expressed (1). 
Vipendutum.center of mass) =~ Lcos a(a@) £- Lsin a(a) ¥ (1) 


Verm = 76 (2) 
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The (1) and (2) can be solved, so axes velocity is utilized (3). 


XS 
ll 


rÅ — Lcosa (a) (3) 


V, = —Lsina (a) (4) 


The Euler-Lagrange formulation has been used to obtain the system dynamic equations such as: 
1) Potential energy: the system gravity energy is presented (5). 


V = P.E(pendulum) = mgh = mgh Lcosa (5) 
2) Kinetic energy: RIP system kinetic energies present (6). 
T = K. Epub + K. Epyx + K. Epy + K. Ependutum (6) 
The RIP rod’s moment of inertia about its COM was presented subsequently as: 
1 
Jem = (=) MR? (7) 
1 1 1 
Im = (EME = (5) ME) = (Yme ® 
After substituting the (3), (4), (7), and (8) in (6), so system kinetic can be obtained (9). 
1 : 1 : ; 1 ; $ 1 ; 
T= (=) Jeq 0? + È) m(r@-Lcosa (&))? + (=) m(—Lsin a(a@ ))? + (=) Jma (9) 
The final Euler-Lagrange formulation can be expressed (10). 
L =T-V = (+) Jeq 6? + (2) mL2a? - mLr cos a(4)(6) + mr20? — mgLcos a (10) 


When 0 and a are the generalized coordinates, so: 


6 (ôL ôL 

= (=) ~ 5a Toutput = Beq 0 (11) 
ô (êL ôL 

Cr er i) 


After linearizing the initial conditions, equations are modified as: 
(Jeq + mr?)6 — mLr = Toutpur — Beg? (13) 
- ml? & — mLr6 — mgLa = 0 (14) 
Then, the motor torque becomes: 


Nm Ng Kt Kg Vm- Kg Km9) (15) 


Rm 


Toutput F 


RIP system state space can be illustrated (16). 


odn i 0 
Ol jo o Oe 0 

aļ_ bd —cG a Nm Ng Kt K 

ale ee Olole (o) 
äl jo © = oflal | omy“ 


2 
Here, a = Jeq + mr?, b = mLr, c = 4/3 ml?,d = mgL, E = ac — b? and G = =" °F. The system 


specification can be listed as shown in Table 2. 
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Table 2. Typical SRV02 system specification 


Symbol Description Value 
K Motor torque constant (N.m\A) 0.00767 
Km Back EMF constant (V.s\radian) 00767 
Rm Armature resistance (Q) 2.6 
Kg SRV02 system gear ratio (motor — load) 14 (14x1) 
Nm Motor efficiency 0.69 
Ng Gearbox efficiency 0.9 
Beg The equivalent viscous damping coefficient (N.m.s\radian) 1.5 e-3 
Jeq Equivalent moment of inertia at the load (Kg.m’) 9.31 e-4 


After substituting the SRV02 system parameters values, then the system state space can be illustrated as: 


Å 0 0 1 07/6 0 
ál _ lo 0 0 ailla 0 
ð -| 39.32 -14.52 o|lòl|* |2554 (17) 
al lo 81.78 -13.98 ollel l24.59 
jooo 0 
Y=[5 4-0 Olle! “lg (18) 


3. CONTROL SYSTEM SCHEME 
3.1. PI-PD and FOPI-FOPD controllers design 

While the typical PID controller has limits in regulating such systems, the PI-PD controller structure 
offers an outstanding four-parameter controller for handling integrated, unsteady, and resonant systems to 
steady state adjustment [19]. PI-PD controller is a mathematical model-initiated theory. High-quality 
fulfillment can be gained at unstable and integrating processes [20], [21]. Figure 2 illustrates the PI-PD 
controller structure as PI and PD parts. Placing the system poles in required places can be achieved by 
interior feedback controller PD, on the other hand, the exterior loop is responsible for controlling the system 
after the interior loop action [22]. So, the PI-PD controller design has features of the traditional PID 
controller. In Figure 2, Gp(s) and D(s) present the controlled system and disturbance [23]. When Gp(s) is: 


Gp(S)= ZO (19) 


PP(s) 
PD and PI can be expressed as Cpp (S) and Cp;(s), sequentially: 


Crp(S) = Ky + Ka(s) (20) 


Ki Kps+ Kj 
S 


Cp (S) = Kp + — = (21) 


Figure 2. Single input single output control system with PI-PD controller [24] 


Where N,,(s) is the numerator and D,(s) is the denominator polynomials of the system transfer 
function, K; is the integer gain, K, is the feedforward gain, Kg is the derivative gain, and Ky is the 
proportional feedback gain. The approach design for FOPI-FOPD controller is the same as the integer PI-PD 
one, except the containing fractional integral and fractional derivative components. In addition to Kp, Kj, Kr, 
and K4, u order of differentiator, A order of integrator represents the additional factors, which maximize the 
controller degree of flexibility to improve performance and achieve more robustness for the controlled 
system [25], [26]. The (22) and (23) illustrate the FOPI and FOPD mathematical description [27]. 
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Crp (S)=Ky + Kas" (22) 
Ki 
Cp,(S) = Kp + a (23) 


A PI-PD controller can be conducted by choosing A = 1, u = 1 [28]. Optimal performance can be 
achieved by fine-tuning the controller’s parameters. GWO technique has been used in this work to acquire the 
best minimizing error rate and reaching excellent stability by requiring the suitable controller parameter [16]. 
The integral time square error (ITSE) performance criteria have been applied as a cost function (24). 


ITSE=J, te(t)*dt (24) 


The two loops of FOPI-FOPD and PI-PD controllers are designed to stabilize the RIP system in a 
SIMO design method illustrated in Figure 3 [29]. So the upright pendulum position and arm angle stabilization 
can simultaneously be achieved [30]. The SIMO structure of the RIP system requires the designing of two 
various controllers [31]. The first controller controls the arm, while the second controls the pendulum [32]. 


FOPI/PD 


Arm + 
Setpoint Rotary 
inverted 
pendulum 
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Pendulum 
Setpoint 


FOPI/PD 


Figure 3. RIP SIMO controlled system illustration 


3.2. Optimization technique 

As control systems cannot accurately predict plant structure in the control system due to parameter 
uncertainty and uncertainty in plant models, soft computing based evolutionary algorithms can be used [33]. 
This work exploits the GWO technique to reduce the ITSE and acquire the best values of FOPI-FOPD and 
PI-PD controller’s factors [34], [35]. GWO is a conceptual program miming the grey wolf hunting strategy 
and social hierarchy. In nature, grey wolves live in a pack and have a stringent social hierarchy. On average, 
5-12 participants in the group [36], [37]. Of specific weal is that they have a very rigorous social dominant 
hierarchy, as shown in Figure 4. The hierarchical social structure consists of four levels of gray wolves. 
The leader, alpha (æ), is at the top of the pyramid. The alpha makes most choices about hunting. Followed by 
beta (£), which supports the alpha in making decisions and carrying out other pack duties. The delta (8) 
wolves provide the information about the wolves’ pack to alpha and beta. Omega (w) is at the bottom of the 
pyramid, which always has to succumb to the other wolves [38], [39]. 


S/N 
fs \ 
A Te 
fs 


Figure 4. Grey wolf hierarchy 


Grey wolves hunting can be presented as a theoretical equation. The alpha, beta, and delta have 
more knowledge than others about geography. So, we preserve the top three possibilities we have found so 
far and force the remaining search factors to update their placings in consideration of the best search factor’s 
position [40], which can be shown: 
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Dy = aX -4.o -Ae - 2 (25) 
X; = Xa — Ay.(Da).X2 = Xp — Az.(Dg), X3 = Xs — A3-(D5) (26) 
Ž(t +1) = 4a (27) 


3 


Where Xi, Xz, X3 are the best solution of (@), (£), and (ô), respectively. Ay, Az, Az, Cr, C,,and Cz are the 
coefficient vectors, respectively, which can be presented as (A= = 2d.r, and C= 2. r2). Where, rı and r, is a 
random numbers between (0-1) and & defined as linearly decreasing coefficient is changed from (0-2) [41]. 


4. RESULTS AND DISCUSSION 

The fractional and integer controllers have been applied to the RIP system to control the pendulum rod in 
the vertical position. The Simulink block diagram of the controlled system in Figure 5 shows the design of the 
controllers under the SIMO approach and a comparison between their results. Figure 6(a) illustrates the pendulum 
response and Figure 6(b) shows the unit step response of the arm for both controllers. Figure 6(c) illustrates the 
control signal behavior for both controllers. To achieve a fair comparison, the same parameters are used for both 
controllers with additional parameters of the fractional controller (A and p), which improve the system 
performance. It is obvious that the performance of the FOPI-FOPD controller is better than the PI-PD controller. 

From Figure 6, the arm enters the state of equilibrium in less than one second without an overshot 
for the FOPI-FOPD controller, while the PI-PD controller takes about two seconds to achieve the set point 
with a 13% overshoot. So, the pendulum rod entered the balancing state when the FOPI-FOPD controller was 
applied faster than it entered the equilibrium state when the PI-PD controller was applied. Although the 
swinging up of FOPI-FOPD is more than PI-PD, it did not exceed 0.04 rad. Figure 7(a) and Figure 7(b) 
illustrate the controller’s response behavior to the disturbance. As a result, the FOPI-FOPD controller was 
more robust than the PI-PD controller in the disturbance test. 

The results of using the FOPI-FOPD and PI-PD controllers for RIP were illustrated in the preceding 
paragraphs. According to these findings, the time response characteristics acquired by the FOPI-FOPD 
controller are more suitable than those achieved by the PI-PD controller. The controller’s outputs are listed in 
Table 3. It indicates that the FOPI-FOPD controller can acquire a low control effort in addition to achieving 
desirable time response parameters, in terms of rise time (tr), settling time (ts), and overshoot ratio (Mp). 


Arm setpoint Disturbance 


Pendulum Angle (rad.) 
Control Signal (volt) 


Ri 
FOPI-PD Control Signal 
----— PI-PD Control Signal 


3 2 3 2 3 4 5 
Time (sec.) Time (sec.) Time (sec.) 


(a) (b) (c) 


Figure 6. Simulation results of the balance control with the FOPI-FOPD and PI-PD: (a) pendulum response, 
(b) arm response, and (c) control signal behavior 
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Pendulum Angle (deg.) 
Arm Angle (deg.) 


0 1 2 3 4 5 0 1 2 3 4 5 
Time (sec.) Time (sec.) 


(a) (b) 


Figure 7. Simulation result after applying disturbance: (a) pendulum response and (b) arm response 


Table 3. Controller’s performances comparison 


Controller Arm Pendulum 
tr (sec.) ts (sec.) Mp(%) tr (sec.) ts (sec.) Mp (%) 
FOPI-FOPD 0.422 0.781 - 0.172 1.174 46 
PI-PD 0.591 2.153 13 0.271 1.710 23 


A comparison of the proposed FOPI-FOPD controller’s performance with the performances of 
controllers created in earlier works (Table 4 and Table 5) shows its efficacy. It indicates that the FOPI-FOPD 
controller achieves the best transient response characteristics among previously studied controllers. The first 
comparison with model-free backstepping (MFBS) controller [42]. The same initial condition is considered 
for the FOPI-FOPD controlled system, as shown in Figure 8(a) and Figure 8(b). 

A second comparison was presented with full-state (FSF) controller [13]. The same set point of the 
arm is taken and applied to the FOPI-FOPD controlled system, as shown in Figure 9(a) and Figure 9(b). 
The results proved that the proposed controller is more efficient than the FSF controller. 


50 


Arm Angle (deg.) 
Pendulum Angle (deg.) 


---- Ref 
FOPI-PD 


oO 1 2 3 4 5 0 1 2 3 4 5 
Time (sec.) Time (sec.) 


(a) (b) 


Figure 8. Simulation result after applying the same initial condition of search [42]: (a) arm response and 
(b) pendulum response 


Table 4. Stabilization performance with previous work comparison 


Controller Arm Pendulum 
tr (sec.) ts (sec.) Mp(%) tr(sec.) ts (sec.) Mp (%) 
FOPI-FOPD 0.266 0.403 — 0.059 0.713 23 
MFBS [42] 0.291 0.513 1 0.161 0.824 23 


Table 5. Proposed controller versus previous work tracking performance comparison 


Arm Pendulum 
Controller tr (sec.) ts (sec.) Mp (%) tr(sec.) Deviation (degree) 
FOPI-FOPD 0.442 0.781 - 0.172 -1.7-0.75 
FSF [13] 0.925 1.107 = 0.181 -2 —0.82 
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Figure 9. simulation result after applying the same setpoint of arms set of search [13]: (a) arm response and 


5. 


The 


(b) pendulum response 


CONCLUSION 
The SIMO designs of FOPI-FOPD and PI-PD controllers for the RIP system are given in this paper. 
GWO approach was used to determine the optimal parameters for both controller gains. The time 


response requirements and ITSE restrictions were incorporated into the suggested cost function to satisfy the 
needs of robustness and time response characteristics. The FOPI-FOPD was designed to improve system 
tracking. The suggested robust FOPI-FOPD controller was proved capable of driving the system, and the 
GWO approach was used to make the controller computationally efficient. Furthermore, the time response 
parameters achieved by the FOPI-FOPD controller are preferable to those acquired by the PI-PD controller 
and controllers developed previously. 
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